4.7.3 Flow with pressure gradient The approximate method based on von 
Karman’s momentum integral equation was extend€d by Pohlhausen to cover 
the case of two-dimensional laminar flow with pressure gradient. We choose a 
coordinate system in which x denotes the distance measured along the surface 


(see Fig. 4.20) and y denotes the normal distance from the surface. The 
momentum integral equation (4.103) is written in the form 


U — + (26+ 6%) U — =—.- (4.113) 

dx dxi p 
In writing down an expression for the velocity distribution, it can no longer 
be assumed that the velocity profiles at varying distances downstream of the 
leading edge are similar. In particular, since dp/dx # 0, we must allow for the 
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existence of velocity profiles with and without a point of inflection. Therefore, 
we assume that the velocity distribution may be expressed as a fourth-degree 
polynomial equation in 7 (= y/5), and we write 


ee 2 3 4 
m aii oto elke EN +A3n~+Agn (4.114) 


where the coefficients may vary with x. To evaluate the five coefficients, we 
stipulate the five boundary conditions 


ðu ldp UdU 


= (0) : u=0, — = — — = — — Á— 
= dy? pdx v dx 
r o ðu 3u 
= u = 3 — m Å——— = 
ý dy ay? 


The first boundary condition at the surface obviously requires that Ag = 0. The 
second boundary condition at the surface ensures that the velocity profile 
exhibits a point of inflection when dp/dx > 0 but has no inflection when 


dp/dx <0. 
We now define the dimensionless parameter 
ô? dU 
A= —— (4.115) 
v dx 
and the four remaining coefficients in Eq. (4.114) are then found to be 
A,=2+ A à 
64, R 
A 
A,=—2+—>, Ag=l1-——- 
> j ai 6 


Therefore, the velocity distribution may be written as 
u A 
Wael ae (m — 31? + 3n? — 74) 
= F(n) + AG(n) (4.116) 
where 
F(n) = 27 — 27? + n* = 1 — (1 +7) (1 — 7)? 
and 2 (4.117) 
Gn) = $ — 39? + 3n? — 74) = paths 7)’. 
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The velocity distribution expressed by Eq. (4.116) is seen to be a 
combination of two distributions. The first component, F(n), is the polynomial 
distribution derived in Example (4.4) for zero pressure gradient. The second 
component, G(n), represents the distortion of the velocity profile associated 
with the pressure gradient, the degree of distortion depending on the magnitude 
of A. Thus, A may be regarded as a shape factor. The resulting velocity 
distributions for various values of A are plotted in Fig. 4.21. It follows from the 


Cle 


Fig. 4.21 Effect of pressure gradient on velocity distribution 


definition of A (Eq. (4.115)) and the fact that dp/dx = —pU(dU/dx) that the 
pressure gradient and the shape factor are related by 


A=-——- (4.118) 


The calculated velocity distributions indicate that the condition (du/dy)g = 0 
(ie., separation) is achieved when A = —12. For A > +12, the calculation 
predicts values of u/U > 1, which have no physical significance in steady flow. 
Consequently, we restrict our attention to the range —12 <A < +12. 

We can now obtain 5*, 6, and 79 in terms of ô and A only. Thus, 


1 1 
ô* u 
=f (: - =) dyn = f [1 — F(n) — AG(n)] dy 
0 0 
nee (4.119) 
“10 120 i 
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1 
9 
ron f [F + AGIU — F) — AG()] dy 
o 


37 A A2 (4.120) 
~ 315 945 9072 ; 


sage coat 
ag dy/o è a Jo 


et a (4.121) 


and therefore 


The shape factor A is still an unknown quantity since it contains the unknown 
thickness 5. If we multiply each term in the momentum integral equation 
(4.113) by 6/pU, we obtain the non-dimensional form 


U6 dé &* \ 8? dU 8 
+ (FEE, (4.122) 


We now define 


62 dU 
K= — — (4.123) 
v dx 
and 
82 
yA a (4.124) 
y 
so that 
K zoe. (4.125) 
= re x 


These substitutions and the analysis which follows represent an improvement on 
the original Pohlhausen method and is due to Holstein and Bohlen. It is seen 
from Eq. (4.123) that K is a second shape factor which is related to 0 in 
precisely the same way that A is related to ô. 

From the definitions of A and K, it appears that 


K AV (= A A Ph (4.126) 
Fy 315 945 9072) ` i 
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Writing 


oA 

ôt 8/5 10 120 

6 O86 Saat A AAS =f,(K) (4.127) 
315 945 9072 


799 “of §-(2+2)(- A A2 
6 /\ 315 945 9072 


#)- =f,(K) (4.128) 
and differentiating Eq. (4.124) to give 
aS, e a (4.129) 


substitution of Eqs. (4.124)-(4.129) into the momentum integral equation 
(4.122) yields 


U dz 
Sae +f, (K)) K =f2(K) 


or 
dZ 
Oe = 2f,(K) — 2K[2 + f,(K)] = F(X). (4.130) 


Therefore, the momentum integral equation can be written in the abbreviated 
form 


dZ F(K) 
=a (4.131) 


where K is given by Eq. (4.125). 

Equation (4.131) is a first-order non-linear differential equation in Z. Since 
K, f,(K), f2(K), and F(K) can all be written in terms of A only, they may be 
calculated for various values of A. Some of the calculated results of Holstein and 
Bohlen are summarized in Table 4.2. 

167 


Table 4.2 Functions for calculation of the laminar boundary layer with 
pressure gradient (Holstein and Bohlen) 


4.7.4 Boundary layer calculations The calculation of the boundary layer 
parameters can be conveniently undertaken as a step-by-step analysis using the 
finite-difference form of Eq. (4.131) 
F(K;) 
AZ, j= ——Ax;,; (4.132) 
U; 

this giving the increment of Z for the j-th step of the analysis, based on the 
calculated values of K and U at the end of the preceding step (Fig. 4.22). The 
calculation begins at the stagnation point and proceeds along the surface until a 
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point of separation is reached. The velocity distribution, U(x), outside the 
boundary layer is assumed to have been obtained by a potential flow analysis. 

To commence the calculation, we need to know the initial values of Z and 
dZ/dx at the stagnation point. From Eq. (4.125), we see that 


Ko 0-0770 
(dU/dx)y (dU/dx)o 


but if we attempt to use Eq. (4.132) to determine (dZ/dx) 9, we obtain an 
indeterminate result since both F(Kg) and Up are zero. However, if we take the 


limiting value 
Js 
—] = lim |— 
dx o x+0 ôx 


dZ (d? U/dx?), 
ae, Be PAE ee a 4.134 
(=), 00652 IOS SER 


Having decided on the size of step (Ax) to be used, the calculation may 
proceed with the following routine: 


1. Calculate Zo and (dZ/dx)o using Eqs. (4.133) and (4.134), respectively. 
2. Calculate 


Zz (4.133) 


it is found that 


3. (i) Obtain A from Table 4.2. 
(ii) Calculate 6 from Eq. (4.123). 
(iii) Obtain 6* from Table 4.2 (fi (K) = 5*/8). 
(iv) Obtain To from Table 4.2 (f2(K) = T09 /uU). 
(v) The velocity distribution can be calculated from Eq. (4.116) if required. 
4. The second step commences using Eq. (4.132), 


F(K;) 
ESEN 
152 U, 1,2 
Then 
Zz =Z; + AZ; 2 
and 
dU 
K,=Z, — 
2 
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5. The values of A, 9, 5*, Tọ, and the velocity distribution are obtained as 


before. 
6. The third step commences with the use of Eq. (4.132), 


F(K2) 
Uz 
2 Che. 


AZ, 3= Ax 3 


The accuracy of the calculation naturally depends upon the number of steps 
chosen; this choice may well be dependent on the availability of computer 
facilities. Adequate calculations can, however, be performed with the aid of a 
desk calculating machine. The method leads to accurate predictions in regions of 
negative pressure gradient. With dp/dx > 0, accuracy deteriorates as the 
magnitude of the pressure gradient increases. 


Example 4.5 Using a step-by-step calculation, determine the distribution of 6, 
5*, and 7ọ on the surface of a circular cylinder and plot the boundary layer 
velocity distributions. 


The potential flow solution for this case was considered in §3.3.3, which gave 


U=2U, sing. 
Then 

dU 1dU 2 

de R&B RS 
and 


d7U 1 d? U 2 
de? Rag? RP 
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(v) The velocity distribution is calculated from Eq. (4.116). 
4. For the second step: 


(K,) 
Z1,2= RAd¢, 2 
0- 5 R 
„ TA Z 0-021 
2U,,sin 5° 57-3 2U, 
Z,=Z,+AZ 0-0791 sé 
2=Z, 1,2= 3g 


dU 
K, =Z, 2 = 0-0791 cos 10° 
2 


= 0-0767. 


5. The values of A, 8,8*, and 79 are then obtained as for the first step. 

6. The procedure is continued until a point is reached for which 79 = 0 (i.e., 
separation). The results of the complete calculation are plotted in Fig. 4.23. 
The results of the exact solution are also included for comparison. 
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————Exact solution 


Fig. 4.23 Approximate calculation of the laminar boundary layer on a circular 
cylinder 


